Amplitudes of Thermal and Kinetic Sunyaev-Zel'dovich Signals from Small-Scale 

CMB Anisotropics 
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While the arcminute-scale Cosmic Microwave Background (CMB) anisotropics are due to sec- 
ondary effects, point sources dominate the total anisotropy power spectrum. At high frequencies 
the point sources are primarily in the form of dusty, star-forming galaxies. Both Herschel and 
Planck have recently measured the anisotropy power spectrum of cosmic infrared background (CIB) 
generated by dusty, star-forming galaxies from degree to sub-arcminute angular scales, including 
the non-linear clustering of these galaxies at multipoles of 3000 to 6000 relevant to CMB secondary 
anisotropy studies. We scale the CIB angular power spectra to CMB frequencies and interpret 
the combined WMAP-7 year and arcminute-scale Atacama Cosmology Telescope (ACT) and South 
Pole Telescope (SPT) CMB power spectra measurements to constrain the Sunyaev-Zel'dovich (SZ) 
effects. Allowing the CIB clustering amplitude to vary, we constrain the amplitudes of thermal and 
kinetic SZ power spectra at 150 GHz. 



I. INTRODUCTION 

The Cosmic Microwave Background (CMB) radiation 
represents a fundamental observable for Cosmology and 
at present the most powerful one for the investigation 
of several open questions, such as the nature of inflation 
or primordial non-gaussianity. In two years the Planck 
satellite [l| will provide a measure of the anisotropics 
of the CMB with a precision never reached before, that 
will allow a highly precise determination of the standard 
cosmological parameters as well as major constraints on 
some non-standard Physics processes. 

The observation of CMB anisotropics is affected by 
several systematics and secondary effects due to the fact 
that the CMB is not the only source of emission in the 
microwave frequencies and to the formation of structures 
between the observer and the last scattering surface [ll^ . 
The great accuracy of future data requires a compelling 
description of these effects, in order to separate the differ- 
ent contributions to the anisotropics and to distinguish 
primordial and secondary effects (see 0). The observ- 
able used to extract most of the cosmological information 
from the CMB is the angular power spectrum Ct- Sec- 
ondary effects or unresolved foregrounds provide a con- 
tribution to the observed Cf. In order to obtain an unbi- 
ased determination of the cosmological parameters from 
CMB maps it is necessary to correctly describe possible 
contaminations. On the other hand, both contaminants 
and secondary effects themselves contain certain cosmo- 
logical and astrophysical information, especially on the 
formation of structure at late-times and the large-scale 
structure of the Universe, so that the separation of these 
components from primordial CMB fluctuations becomes 
an important science goal on its own. The Galactic emis- 
sion and radio point sources are typical examples of fore- 
ground contamination in CMB maps. While the bright 



sources detected in maps can be removed with a suit- 
able mask before the estimation of the angular power 
spectrum, unresolved point sources will contribute to the 
total anisotropy power spectrum Cg. 

The Sunyaev-Zel'dovich (SZ) effect 0, caused by the 
Compton scattering of the CMB photons by the electrons 
in the Universe, is a well-known secondary anisotropy 
studied by a variety of experiments. The SZ effect con- 
tains cosmological information, since the angular power 
spectrum of secondary temperature anisotropy arising 
from SZ scattering depends on both the gas distribution 
in galaxy clusters and on the amplitude of the matter 
density fluctuations trg j8l-[ll|. Since the SZ thermal effect 
has a unique spectral signature relative to the CMB ther- 
mal spectrum, the SZ signal can be distinguished from 
primary CMB anisotropics and other foregrounds using 
observations at multiple frequencies across the SZ null 
at ^ 217 GHz [l^]. Such a separation, however, is not 
feasible with kinetic SZ effect associated with peculiar 
motions of the electrons scattering the CMB Q as the 
signal has the same spectrum as that of the CMB. Even 
for the SZ thermal effect, in realistic experiments, the 
main obstacles that limit a clear detection of the thermal 
signal comes from uncertainties in the modeling of the 
kinetic contribution and the difficulty of separating SZ 
effects from clustered point sources. 

The contribution of unresolved point sources and SZ 
effect is best seen on small angular scales where they 
dominate the total CMB angular power spectrum. The 
use of data at these small angular scales is hence becom- 
ing decisive in the analysis of CMB data. In this work we 
analyze large (. data (up to I ^ 9000) from the South Pole 
Telescope (SPT) at 150 and 220 GHzJlland from At- 
acama Cosmology Telescope (ACT) [iJ: 150 GHz 
combined with Wilkinson Microwave Anisot ropy Probe 
(WMAP) data after 7 years of observation [1^1 to put 
constraints on the two SZ effects. 
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In previous studies significant limitations came from 
uncertainties associated witli clustering of dusty star- 
forming galaxies (DSFG) that contribute to high- 
frequency CMB data. Such clustering has now been mea- 
sured with both Herschel p7| and Planck [l| experiments. 
In the context of CMB studies, Herschel measurements 
are most useful as they probe the DSFG clustering down 
to sub-arcminute angular scales at scales well matched to 
arcminute scale CMB experiments while Planck measure- 
ments are limited to scales greater than 5' or ^ < 2000. In 
this work we describe the clustering of unresolved point 
sources we used the same template of [13 , where the au- 
thors reported a detection of both the linear clustering 
and the non-linear clustering at a few arcminute scales, 
corresponding to £ ^ 4000. 

We perform a Monte Carlo Markov Chain (MCMC) 
analysis constraining both the thermal and the kinetic 
terms of the SZ effect together with the the Poisson and 
clustering corrections due to unresolved point sources, 
including radio sources at lower frequencies such as 150 
GHz. In the next Section we recall more details on 
the contribution of the SZ effect and foregrounds to the 
CMB anisotropy power spectra. Section 3 describes the 
parametrization and the templates we used to model the 
foreground and the SZ contamination to the CMB angu- 
lar anisotropy power spectrum. In Section 4 we show the 
results and conclude with a summary. 

II. PARAMETRIZING SZ EFFECT AND 
FOREGROUNDS 

In this Section we briefly describe the adopted 
parametrizations and templates for the Sunyaev- 
Zel'dovich effect, unresolved extragalactic point source 
foregrounds and lensing. 

Sunyaev-Zel'dovich thermal and kinetic effect. 

The SZ effect has two different contributions, one from 
the thermal motion of the electrons (thermal SZ effect - 
tSZ) and one from the bulk motion of the electrons rela- 
tive to the CMB (kinetic SZ - kSZ). The former contribu- 
tion has a distinct frequency dependence, while the kSZ 
effect causes only a Doppler shift of the incident CMB 
spectrum retaining the black body shape. The total SZ 
signal in a generic direction h is then given by (see for 
example Section 2 in [l3|): 

/(a^o) 

with X = hiz/ksTcMB and f{x) = 2 — x/2 tanh(2;/2). 
Here ATq^^ is the expected thermal contribution at fre- 
quency vq. From f{x) it can easily seen that the ther- 
mal SZ effect vanishes at ~ 218 GHz. We model the SZ 
contributions to the anisotropy angular power spectrum, 
relative to a template power spectrum, as 

nSZ,lJ A /(^l) fi'^j) ntSZ I A nkSZ (ryS 

= ASZTT— T77— T-t^O/ + ^kSzi^O.f : i^J 



where Di = i{£+l)Ci/2Tr and Dq ^ is the template spec- 
trum for either thermal or kinetic SZ. In this work we 
consider the SZ templates from [isj . computed by trac- 
ing through a dark matter simulation and processed to 
include gas in dark matter halos and in the filamen- 
tary intergalactic medium. The thermal SZ template 
describes the power from tSZ temperature fluctuations 
from all clusters for a Universe normalized with am- 
plitude of matter fluctuations as ~ 0.8. In particu- 
lar we use the 'standard' model of [ll|, that was first 
described in [l^, and assuming a ag scaling given by 
Dg^f oc (ct8/0.8)^'^ as found in [l^. For these templates 
the reference values a.t £ = 3000 are D^^^^^qq^ ~ 8.9iiK^ 
and i^o'fcsooo - "^-^l^K^ 

Foregrounds from unresolved extragalactic point 
sources. 

The foregrounds contribution to the CMB power spec- 
trum at arcminute angular scales arises essentially from 
unresolved extragalactic point sources. These sources 
provide two contributions, a Poisson term due to the ran- 
dom discrete distribution and a clustering term account- 
ing for the large-scale distribution of the sources. We 
assume the Poisson term as constant in Ce, modeling it 
as Df°'^^ = ^Poissi'o.f where D^°i^^ = (£/3000)2. The 
clustered term can be similarly expressed as — 
^ciust(i^i, where -Dq'^'^' is the point sources clus- 

tering template and ^ciust(t'i, ^j) encodes the frequency 
scaling (see section 3 for further details) . Contribution to 
point sources comes from radio point sources and dusty 
star-forming galaxies (DSFG). At 220 GHz the main 
point source contribution is mainly DSFGs while at 150 
GHz the point sources arc primarily radio sources with 
a synchrotron spectrum. We therefore neglect the clus- 
tering of radio sources and assume that the contribution 
from radio sources is essentially described only by a Pois- 
son behavior. For clustered DSFGs we adopt the tem- 
plate from [TtI where the authors reported a detection 
of both the linear clustering and the excess of cluster- 
ing associated with the 1-halo term at arcminute scales. 
Those data are from the Herschel Multi-tiered Extra- 
galactic survey (HerMES) [2^, taken with the Spectral 
and Photometric Imaging Receiver (SPIRE) onboard the 
Herschel Space Observatory [2l| . 

CMB Lensing 

It is well known that gravitational lensing of CMB 
anisotropics by large-scale structure tends to increase the 
power at small angular scales (see [22| for a recent re- 
view). A proper calculation of this effect is hence neces- 
sary in order to prevent an incorrect estimate of the fore- 
grounds and SZ parameters. The calculation of lensed 
CMB spectra out to £ = 9000 is prohibitively expensive 
in computational time. Instead, we approximate the im- 
pact of lensing by adding a fixed lensing template D'f"''' 
computed by running cernib [23| with and without the 
lensing option and taking the difference between these 
spectra. In this run the cosmological parameters of the 
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ACDM model arc fixed at the best fit values WMAP7. In 
[l3| it has been estimated that the error due to this ap- 
proximation is less than 0.5^K^ at ^ > 3000 and is hence 
negligible with respect to secondary and foregrounds con- 
tributions. The lensing contribution is clearly frequency 
independent. 



impact of this assumption. Finally we consider a sin- 
gle clustered point sources parameter, ^ciust, scaling the 
contribution at different frequencies using the relation of 



III. ANALYSIS METHOD AND DATA 
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We place constraints on the cosmological parameters 
and on the SZ and foregrounds parameters using the 7- 
years WMAP data in combination with the SPT data at 
150 GHz and at 220 GHz, and the AGT data at 148 GHz. 
The SPT and ACT datasets are necessary to analyze the 
smaller scales of the power spectrum where point sources 
and SZ are dominant. For the SPT data we select the 
single frequency 15 x 15 blocks from the full 45 x45 covari- 
ance matrix provided by the SPT collaboration (see [l3|). 
neglecting the correlation between different frequencies. 

We use a 6-parameter flat-ACDM cosmological model 
to describe primary CMB anisotropics and reionization: 
the baryon and dark matter physical energy densities 
Vlbh? , Vlch? , the reionization optical depth r, the ratio 
of the sound horizon to the angular diameter distance 
at the decoupling 0, the amplitude of the curvature per- 
turbation As (with flat prior on logAg) and the spec- 
tral index tIs ; these two last parameters are both defined 
at the pivot scale fco = 0.002 hMpc~^. In addition to 
the standard cosmological parameters we include the SZ 
and foreground parameters described in the previous sec- 
tion. We perform a Monte Carlo Markov Chain analy- 
sis based on the publicly available package cosmomc |24{ 
suitably modified to account for the additional parame- 
ters, with a convergence diagnostic based on the Gelman 
and Rubin statistics. When estimating parameters with 
point sources and SZ included, the total CMB anisotropy 
spectra arc three, one for each frequency plus the cross- 
correlation term, because of the frequency dependence of 
the secondary anisotropics. 

In order to study the stability of our results on 
the assumed paramctrization, wc perform three differ- 
ent analysis, both with 6 additional parameters describ- 
ing SZ effect and foregrounds, but considering different 
paramctrizations. 

First case: "runl". 

In the first case, that we define as "runl" in what 
follows, we consider the SZ effect parameters ^tsz and 
^ksz, the Poissonian contribution ^po^gg and A'^^^^^ and 
the Poisson contribution for the 150x220 GHz cross- 
correlation A-pQjgg . The use of the ^poigg extra parameter 
for the cross correlation of Poisson point source can be 
justified from the possibility that the contribution at one 
single frequency comes from more than one point source 
population and that the two channels are not fully cor- 
related. This possibility has not been considered in pre- 
vious analyses and is therefore important to evaluate the 



with vq = lOOcm^^, following recent results from Planck 
[3 . In what follows we refer to this scaling as " Gispert" 
scaling. 

In summary, the spectra in "runl" are defined as: 
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The thermal SZ effect is negligible at 220 GHz. The 
contribution of the thermal SZ effect to the cross- 
correlated power spectrum may not vanish in presence 
of a spatial correlation between IR sources and the clus- 
ters that cause the thermal SZ. Nevertheless as showed 
in [l^ the effect of this correlation is negligible for the 
SPT data (see par. 7.4 in [11] for further details). We 
hence do not consider this contribution when fitting the 
data. 

Second case: "run2" . 

In the second analysis, to which in what follows we 
refer as "run2", we assume full correlation between the 
Poisson point sources signal at 150 and 220 GHz as done 
in previous analyses, i.e. we fix the cross amplitude of 
Poisson point sources at the square root of the product of 
the ampUtudes at 150 and 220, ^^.^gg - V^Poiss^PoL- 
Moreover, we don't scale the clustered point sources tem- 
plate and we use instead two different parameters for 150 
and 220 GHz. This second analysis is more similar to the 
one presented in [Tsj . however we point out that while 
here we consider the amplitudes at different frequencies 
as free parameters, [l3l | considered the amplitudes at one 
single frequency and one common frequency spectral in- 
dex for clustered and point sources as free parameters. 

In this second case the foreground spectra are defined 
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Third case: "runS" . 

Finally we combine 150 GHz data of SPT and ACT, 
using separate parameters for ACT and SPT both for 
clustered and Poisson point sources, to account for the 
different masking thresholds of the point sources. In this 
case wc have: 
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FIG. 1. Joint two-dimensional posterior probability con- 
tours showing 68% and 95% C.L. constraints on D^'^^qqq and 
DEsooo from ACT 150 GHz data (red) and SPT all frequen- 
cies data (blue) for the run2 case. 



IV. RESULTS 

In Tabic |T] wc report the mean values of the cosmo- 
logical parameters and their 68% C.L. uncertainty from 
SPT data at 150 and 220 GHz for the "runl" and "run2" 
analyses, while in Table [H] we list the mean values of 
the cosmological parameters and their 68% C.L. uncer- 
tainty from ACT data ("run3") combined with SPT data 
at 150 GHz. In order to facilitate the comparison with 
other works present in the literature we also translate 
the constraints on the foregrounds amplitudes in to the 
foreground power spectrum at £ = 3000, I?£=3ooo- Since 
a significant correlation exists between thermal and ki- 
netic SZ and since the kinetic SZ is predicted to be small, 
we also perform an analysis by fixing D'^^qqq = 2 i^lK^ . 

We find that for the "runl" case the thermal SZ 
anisotropy amplitude is -D^lsooo ^ 2.2 ± 1.5 fiK^ . 
While a Icr indication for SZ is present our result 
is less significant than the one reported by [l^l with 
-^Elooo = 3.2 ± 1.3 fiK^ i.e. with a thermal SZ de- 
tection at at more than two standard deviations. The 
result on the kinetic SZ component are compatible, with 
-D^liooo ^ 2.7 ± 1.9 /lA'^ at 68% c.l. from our analysis to 
be compared with i'fffooo = 2.4 ± 2.0 /lA'^ from [l!]. 

Although the point sources and SZ parameters do 
not show significant degeneracies with cosmological pa- 
rameters (see also H^), a strong correlation exists be- 
tween Atsz and ^ksz and to a smaller extent between 
^tsz,ksz and Adust- This can be seen in Figure [T] where 
wc show the 2 — D likelihood constraints in the plane 



-C'fciooo ~ -D/'£looo for the "run2" and "runS" case. Fix- 
ing the kSZ term sligthly improves the detection for the 
thermal SZ with i^Elooo = 2-3 ± 1.4 ^A'^ in "runl" but 
still with less significance than the one in [l3| where a 
value of ffilooo = 3.5 ± 1.0 fiK'^ is reported. 

Based on the degeneracy direction of Figure [TJ we 
constrain the sum of the SZ effects at £ = 3000 to be 
-CtFz°°° + 0.5Di^^° = 3.5 ± 1.8 fjK^ to be compared 
4.5 ± 1.0 ^iK^ of [III . These amplitudes are consistent 
but, again, the significance of the detection is worse than 
who found this sum to be 4.5 ±1.0 /iK^. 

The small discrepancy with the results presented in 
[isj comes essentially from the different paramctriza- 
tion used. Adopting a more similar parametrization 
as in the case of "run2" we found Dg^^ooo — 2.9 ± 
2.0 fiK^, i:'^£fooo = 5.5±3.0AiA:2 at 68% c.l., 0^^^°° + 
0-5^ksz°°° = 5-6 ± 2.6 fiK^, yielding a detection for the 
thermal SZ with higher significance. In case of fixed kSZ 
we obtain D^^qqq = 4.7 ± 2.0 fJ.K^, again a more signifi- 
cant detection in better agreement with jT3| . 

The different assumptions in the frequency scaling of 
the clustered point sources component in "runl" and 
"run2" is the main explanation for the difference in the 
results. In "runl" , taking into account Gispcrt scal- 
ing [ill, we find Df^llll" = 39.11 ± 6.79 fiK^, while 
in "run2" , when the amplitude of the clustering point 
sources is allowed to vary, we have Df^^^^ = 47.33 ± 
9.78 ^A'^, that is more consistent with the correspond- 
ing value of Df^llll" = 57 ± 9 reported in [H. A small 
tension therefore exists between the Gispert scaling and 
the data at 220 GHz, resulting also in a worse detcrmi- 
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TABLE I. Mean values and 68% error bars from SPT data at 150 and 220 GHz. Runl case is with only one DSFG clustering 
amplitude allowed to vary and frequency scaling fixed from [2^, consistent with Planck [l|. Run2 case is with two DSFG 
clustering amplitudes allowed to vary, without frequency scaling. 
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TABLE II. Mean values and 68% error bars from ACT data combined with SPT data at 150 GHz 
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nation of the thermal SZ signaL The use of a different 
parametrization of the point source (just amphtudes in 
our case while [13] varies one amplitude and one spectral 
index per component) can explain the remaining differ- 
ences. 

Concerning the Poisson point sources component at 
150 GHz, our results are different when directly com- 
pared to those from [31, both in "runl" and "run2". At 
150 GHz we find Df^^^l^° = 10.03 ± 0.67 ^iK^ in "runl" 
and D^°^° = 10.38 ± 0.63 ^iK^ in "run2", while the 
value in |l| is -DF=30oo° = 7.4 ± 0.6 /iif 2. This difference 
is explained in straightforward terms if we take in account 
that in [l^ radio galaxies are included in their "baseline 
model" with an amplitude ^^^^3000 = 1-28 /iX^ with a 
15% uncertainty. Clustering of radio galaxies is negligi- 
ble, so this radio galaxies term is a Poisson like term of 
the form oc Adding this component, our Poisson am- 
plitudes are consistent with those reported by [l3| within 
Icr. In "runl" at 150 GHz we find °ijodo° = 10-03 ± 
0.67 /i/v while in [l^l the correspondent total Poisson 
contribution at ^ = 3000 is about (8.68 ± 0.69) fiK'^. At 
220 GHz we find Df°}^^^^° = 79.5 ± 4.8 fiK^ in "runl" 
and £>f°'^° = 77.89 ± 4.49 ^i^^ "run2", while the 
value in gsf is I^F^sooo" = 71 ± 5 fiK^. 

We can therefore conclude that the current results pre- 
sented in the literature on the amplitude of the secondary 
anisotropics should be considered with great care since 
there is a clear dependence on the parametrization used, 
on the frequency scaling adopted and on the assumed 
templates. We stress that, a part for small discrepancies 
imputable to differences in the parameterization, all our 
results for the SZ amplitudes from the analysis of SPT 
data both for our "runl" and "run2" cases, are substan- 
tially consistent with the analysis of the same data made 
by [l3|, even if we are finding less tight constraints. Our 
results hence compare in the same way to the recent pre- 
dictions of tSZ power made by the models of [13], [3 
and [2^ confirming that these models overestimate the 
power of the tSZ signal, as already found in [isl ]. 

In Figurc[2]we show the recent CIB power spectra data 
of the Planck collaboration [l[ at 217 GHz and small an- 
gular scale CMB power spectrum data from SPT, at 220 
GHz, with a comparison to scaled measurements from 
[l3] . The Herschel model is shown in terms of the 1-halo 
and 2-halo contributions to the total power spectrum. 
For reference, we also show the model used by [23 at 220 
GHz to describe the clustering of DSFGs, which overesti- 
mated the power at tens of arcminute angular scales and 
above relative to Herschel and Planck DSFG clustering 
measurements. [291 ] used a linear model to analyze their 
data. At small angular scales non-linear effects are not 
negligible and using a linear model to interpret the data 
may lead to a wrong determination of the bias and hence 
to an overestimation of the power at larger angular scales 
(see also discussion in [l|). Instead our model shows a 
good fit of both Planck and SPT CIB. 

In Figure [3] we show the best fit models for each com- 
ponent compared with the SPT and WMAP7 data. In 
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FIG. 2. Planck DSFG clustering data (red points) at 217 GHz 
and SPT data at 220 GHz (white points) compared with the 
combination (solid black line) of the Poisson term (green line) 
and clustering term (red) of unresolved point sources by scal- 
ing the best-fit model to measurements made with Herschel 
at 350 /xm to 217 GHz using the frequency scaling of [2^. We 
show the 1-halo (pink) and 2-halo (orange) contributions to 
the clustering term following [T^]. The dashed lines are the 
220 GHz SPT DSFG power spectrum components from 
which resulted in an overestimate of Planck DSFG clustering 
at ^ < 3000. 



the "runl" , when only one amplitude of clustered DS- 
FGs is allowed to vary when fitting the all frequencies 
SPT data combined with WMAP7 data, we find that 
^ciust = 1-05 ± 0.19. This suggests that the combi- 
nation of [13] model and the frequency scaling for the 
mean CIB is a good fit of the DSFG clustering at lower 
CMB frequencies. Higher precision CMB power spectra 
at 150, 220 and 350 GHz and a direct cross-correlation 
of _ffersc/ie/-SPIRE maps against the CMB will be neces- 
sary to study if fluctuations scale with frequency as the 
mean CIB intensity and to improve overall constraints 
on secondary anisotropics. 

V. CONCLUSIONS 

In conclusion, we provided a new analysis of the fore- 
ground contribution to the CMB data making use of the 
latest ACT and SPT results. Our work is complemen- 
tary to those presented by the SPT and ACT experimen- 
tal teams since we use a different parametrization of the 
foregrounds contribution and different templates. The 
foreground contribution from Poisson point sources at 
220 and 150 GHz is detected with very high significance 
(at more than 15 standard deviations) with no par- 
ticular dependence on the parametrization used ("runl" 
and "run2" cases are giving very consistent results). The 
contribution from clustered point sources is also well de- 
tected at 220 GHz. We have found that current CMB 
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FIG. 3. Contribution to the angular anisotropy power spectrum from point sources and from SZ effect for the best fit model 
of the WMAP7+SPT analysis. Left panel is 150 GHz, middle 220 GIfz and right panel shows the cross spectra. kSZ term 
is the orange solid line and the tSZ term at 150 GHz is the purple line. Green lines are the Poisson terms and blue lines are 
the clustering contributions. The black lines are the total best fit power spectra. Black dots are SPT data and red squares 
are WMAP7 data. The bottom panels show the residual relative to the total model, including primordial CMB and best-fit 
secondary anisotropy amplitudes. 



data favours a larger contribution at this frequency than 
the one expected by the Gispcrt frequency sealing onee 
the data is normalized at 150 GHz. The thermal SZ com- 
ponent is detected at a level slightly above the two stan- 
dard deviations. However a different parametrization of 
the components and the assumption of the Gispert seal- 
ing could bring this detection to about one standard devi- 
ation. The correlation with the kinetic SZ term is present 
in the data despite the multi-frequency approach. More 
data at more frequencies are clearly needed to establish 
a strong detection of the SZ term. 

While our constraints do not improve the results in the 
literature, we have made a significant addition to prior 
studies by firmly establishing the power spectrum of DS- 
FGs that dominate the arcminute scale CMB anisotropics 
at 220 GHz and higher frequencies. This comes from the 



recent Herschel results combined with Planck-confirmed 
frequency spectrum for the CIB mean intensity. In fu- 
ture, additional improvements will come from directly 
cross-correlating the CMB maps against high-resolution 
CIB maps from Herschel; for this a Herschel-SPIKE sur- 
vey at the same large areas as CMB surveys will become 
useful 
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